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We review various aspects of magnetoresistance in (quasi- )twodimensional systems subject to 
an in-plane magnetic field. Concentrating on single-particle effects, three mechanisms leading to 
magnetoresistance are discussed: the orbital effect of the magnetic field - due to inter-subband 
mixing - and the sensitivity of this effect to the geometrical symmetry of the system, the interplay 
between spin-orbit coupling and Zeeman splitting, and the influence of the field on spin scattering 
at magnetic impurities. 

I. INTRODUCTION 

Studies of low-dimensional electron systems subject to an in-plane magnetic field have come to the focus of intense 
attention recently. In several experiments the influence of an in-plane magnetic field on two- and one-dimensional 
electrons in semiconductor heterostructures HJ^] as well as in lateral quantum dot devices H has been investigated. 
The use of the fairly unconventional in-plane field geometry aimed at achieving a stronger magnetic field influence 
on the 2D electron spin, thus, compensating the dominance of orbital effects in most of the known semiconductor 
materials. Then, by manipulating the field-induced spin polarization, one may extract information about the ground 
state properties of interacting electrons. At low temperatures - i.e., in the regime where electron-electron interactions 
open the possibility of a metal-insulator transition or the ferromagnetic Stoner instability - interaction effects coexist 
with single-particle interference effects. The aim of this article is to provide a theoretical overview of possible influences 
of an in-plane magnetic field on single-particle quantum transport phenomena in semiconductor heterostructures, 
quantum wells, and lateral dots. 

Quantum transport phenomena have been the subject of extensive experimental and theoretical investigations 
during recent decades |7|. There exists a broad variety of electronic systems, where spectacular effects produced 
by the interference of electron waves have been observed. These include weak localization (WL) M, leading to the 
low-temperature magnetoresistance in low-dimensional electron structures, such as two-dimensional electron and hole 
gases in semiconductor heterostructures and quantum wells, as well as universal conductance fluctuations (UCF) ]9|-p^] 
in quantum wires and dots. 
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FIG. 1. Back-scattering from an impurity: Contributions due to multiple scattering events. 

Weak localization of electron waves is the result of an enhancement of back-scattering from a scattering center 
in a disordered environment. This enhancement of back-scattering suppresses electrical conductivity of a quantum 
conductor, as compared to the value expected on the basis of a classical Drude formula. The Drude result is determined 
by the mean free path corresponding to a single-impurity scattering cross section whereas multiple scattering events 
as depicted in Fig. [j] are responsible for the enhanced back-scattering. It is generated by the constructive interference 
between pairs of electron waves that scatter from the same impurities, i.e., two electron waves arrive together at a 
particular impurity, scatter from it to follow a random-walk path (in the form of a closed loop) visiting the same 
set of surrounding impurities - cither in clockwise or anti-clockwise direction - and finally scatter into the exactly 
backward direction. An illustration of electron paths involved in such a non-local scattering process is shown in Fig. |^. 
Obviously, the observation of such an interference effect requires phase-coherent propagation along the closed-loop 
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part of the geometrical trajectory. Furthermore, the magnitude of the observable contribution to the conductivity is 
sensitive to the fundamental symmetries of the system, in particular, to the presence or absence of time-reversal (T) 
symmetry. 

The effect of enhanced back-scattering off an impurity inside a disordered metal can be discussed in the same terms 
as the enhancement of back-scattering of waves from a diffusive medium. The probability u>k,-k of a wave to be 
reflected back from a medium after multiple scattering inside is determined as the modulus of the sum of amplitudes 
Aj corresponding to particular paths j, 



Wk-k 



J2[Af(K-k) + Af(k,-k)] 



(1) 



Here, O (O) stands for clockwise (anti-clockwise) propagation along the path. 

In an ideally time-reversible system, the two amplitudes related to the electron propagation along the same ge- 
ometrical loop { f in opposite directions, ^4p(k, — k) and A!p(k, — k), are identical. Thus, AjA®* = |^4p| 2 , though 
for each particular loop Aj contains a possibly large random semiclassical phase of propagation. The cancellation of 
phase factors happens because for an electron described by a T-symmetric Hamiltonian, (a) the phases acquired along 
the same 'free' segment of the path are equal for waves propagating with the wave number k and — k, and (b) the 
amplitudes of intermediate scattering processes, k a — * k& for and — k/, — > — k a for also equal each other at each 
node of the loop j. Therefore, after averaging over different loops which is equivalent to the averaging over disorder, 
one can identify two phase-insensitive contributions to the back-scattering probability, 



= E (K°f + I4T) + 2Re£ (AfAf*) 

3 3 

= w cla * s +6w Wh , (2) 

where w class is the back-scattering probability of a classical particle whereas 5w WL is an additional quantum contri- 
bution due to constructive interference between time-reversed paths. This second term in Eq. (Q) increases the return 
probability, i.e., the probability for an electron to revisit the same scatterer again - which explains the name 'weak 
localization' given to this effect. 




FIG. 2. Enhanced back-scattering caused by the interference of electrons traveling along the same geometrical path in 
opposite directions. 



The length of closed paths contributing to the interference between back-scattered waves is limited by the phase 
coherence time r^,, yielding a maximal length £ max = v-pT^. In practice, this time scale is determined either by 
phase relaxation within the system or the electron escape from a fully phase-coherent conductor to bulk reservoirs. 
Estimates for the weak localization corrections to the conductance read 

[ ln V 2D ' 

Sg WL ~ - ^ x J _L Q1D wire, (3) 
[ 1 OD dot. 

For the case of a dot, by g we understand its full conductance. Furthermore, r is the elastic mean free path, and 
D = Vpr/d (up Fermi velocity, d dimensionality of the system) is the diffusion constant. 

The violation of T-invariance lifts the equivalence between the amplitudes Aj and A'j' and, thus, suppresses 
contributions to the enhanced back-scattering probability from longer paths. An obvious reason for t — > — t symmetry 
breaking is an external magnetic field: thus, provided the orbital electron motion is coupled to the field, it affects the 
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interference effects in the phase-coherent electron transport. For two-dimensional electrons in a heterostructure subject 
to a magnetic held perpendicular to the 2D plane, the difference between the amplitudes A? and A9 arises due to the 
Aharonov-Bohm phases accumulated along the path. For clockwise and anti-clockwise propagation around the oriented 
area Sj, the Aharonov-Bohm phases, SifAB = B z Sj/4>o (4>o flux quantum), have opposite signs. As for random-walk 
trajectories the encircled area Sj is random, contributions from paths encircling a magnetic flux larger than typically 
the flux quantum (\B z Sj\ > (f>o) cancel out after disorder averaging: (A^A'^*) = ( |A^| 2 x exp[2i<5<^ B ] ) — * 0. 

When the orientation of the applied magnetic field lies within the plane of the heterostructure where an effectively 
two-dimensional electron moves (i.e., B z = 0), an analysis of its influence on the weak localization properties requires 
taking into account more subtle arguments than a straightforward consideration of the Aharonov-Bohm effect. In 
fact, a strictly 2D system is insensitive to the orbital T-breaking effect of an in-plane magnetic field. A possibility to 
couple a 2D orbital motion of the electrons to the in-plane magnetic field and, therefore, to break the time-reversal 
symmetry in a 2D system appears only after taking into account p3-17 the finite extent of electronic wave functions 



in the confinement direction of a quantum well as well as subband mixing by the magnetic field |^8|-|lJ and possibly 
by disorder. 



In Chapter III, we discuss this purely orbital effect of an in-plane magnetic field on quantum interference in 2D 
semiconductor structures. As for the confinement of the transverse motion, we also describe the crossover from 
the quantum limit (where the width of the system is of order of the Fermi wavelength) to the semiclassical regime 



corresponding to thin films of disordered and pure metals |22|- [24| , |1 3| , |25|Jl 5[Jl Tfl ■ The calculation aims at identifying 
the time scale tb after which the time-reversal symmetry breaking sufficiently affects the phase of an electron to 
destroy constructive interference. Correspondingly, only paths with lengths shorter than C(B) ~ vftb contribute to 
weak localization. Thus, tb is the time scale that determines the estimates for the weak localization correction to the 
conductivity, i.e., it replaces the phase coherence time in Eq. (|^). Therefore the conductance is field-dependent, and 
the so-called magnetoconductance Ags = 5g wh (B) — Sg WL (0) obtains as 

2 f ln^ 2D, 
h QID wire. 



More details about the multi-subband case can be found in chapter VI 

Another possibility for the in-plane magnetic field to affect the interference of 2D electron waves arises due to the 
coupling between orbital motion and electron spin, which is a pronounced feature of the electronic band structure of 
most III-V semiconductors with the zinc-blend type crystalline lattice and a unit cell without inversion symmetry. In 
a heterostructure or quantum well, spin-orbit coupling for a spin-i carrier is linear, both in the planar momentum 
k and the spin operator a. For a given initial spin-polarized state, it produces an electron spin precession, where 
the precession axis as well as the precession frequency depend on the direction and value of the Fermi momentum 
of the propagating electron. For a diffusing electron, such a precession randomly changes direction, thus, providing 
spin relaxation for an initially spin-polarized particle - the mechanism known as Dyakonov-Perel relaxation f p6| . 
However, spin-orbit coupling alone does not violate time-reversal symmetry, though it does modify the above-described 
interference effects. In particular, strong spin-orbit coupling changes the constructive interference between electrons 
encircling the same closed geometrical path in opposite directions into destructive interference |2l] |29| . To trace the 
origin of this change, let us consider the extreme limit of dominant spin-orbit coupling, where an electron always 
remain in the same chiral state along its orbital motion (for instance, £ = tK- -a = +1). Thus, an electron approaching 
the scattering medium in the initial state |k, £) (which can be viewed as |k, — ►)) ends up in the final backward moving 
state | — k, £) (which can be viewed as |— k, *—))■ This implies a rotation of its spin within the x — y plane by the 
angle \0j\ = (2n + l)ir (where n £ N) while moving along a diffusive path connecting these two states. For two spin-i 
electrons encircling the same closed loop in opposite directions, 9® = —9^ — (2n+ 1)tt. As a result, their interference 
produces a negative contribution to the back-scattering, 



A°A°* = |A°| 2 exp 



9? - ef 



Af\\ (5) 



which increases the quantum conductance as compared to its classical Drude value - the effect known as weak anti- 
localization. 
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FIG. 3. Suppressed back-scattering in the presence of strong spin-orbit coupling. 



Now, the application of an in-plane magnetic field causes a Zeeman splitting of spin-| and spin- J, electron states 



which - in combination with spin-orbit coupling - affects the interference of back-scattered waves. Chapter IV is 
devoted to the description of the interplay between Zeeman splitting and spin-orbit coupling in determining quantum 
transport characteristics (in particular, of lateral dots). 

Interference between different paths not only affects the mean value of the conductance (i.e., (g) — g c \ ass + <5g WL ) 
but also characteristic deviations from that mean value, g = (g) + Sg. Universal conductance fluctuations (UCF) 
represent another interference phenomenon observed in mesoscopic conductors and quantum dots at low temperatures. 
To discuss this effect, let us consider the transmission probability through a disordered piece of metal. For a diffusing 
electron wave, there are many scenarios leading it from the same initial state |k)i e ft on the left to the same final 
state |k') r ight on the right, thus, creating the possibility for interference. Depending on the individual distribution 
of scatterers in the sample, different paths j have different weights - and, thus, the waves following them acquire 
different phases - which makes the transmission probability, 



k,k' 



E^-( k ' k ': 



El^f + ^A*, (6) 



j i 



sample-specific as well as dependent on specific interference conditions, namely the energy (wavelength) of the incident 
electron or the presence of an external magnetic field. In measurements, this manifests itself in a random but 
reproducible dependence of the sample conductance on the Fermi energy E-p and the applied field B z . Furthermore, 
gate voltages (varying the shape of a semiconductor dot) represent an additional control parameter. 
In a phase-coherent system, the variance of conductance fluctuations has a value, 
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(Sg ) = c gcom a x I — j , (7) 

which is universal |Tl||3C|]3l| ] up to a geometry-dependent numerical factor c gcom ~ 1. Furthermore, a reflects the 
fundamental symmetries of the system, expressed by three coefficients, 

The symmetry plays such an important role here as it determines the number of independently fluctuating components 
of the generically random amplitudes Aj related to each particular geometrical path. The higher the symmetry, the 
more correlations are implicit between real and imaginary parts as well as spin components of Aj. In particular, s = 1 
or 2 depending on whether Kramers' degeneracy is lifted or not. When the orbital electron motion is independent of 
the electron spin state, i.e., when the system possesses spin-rotation invariance, the T-symmetric system is described 
by j3 = 1 whereas f3 = 2 if time-reversal symmetry is broken. The value (3 = 4 (with s = 1) represents the case of 
efficient spin-relaxation due to spin-orbit scattering. Finally, an additional symmetry class number E |H| characterizes 
whether a violation of time-reversal symmetry in the spin-sector affects (E = 2 - or not, E = 1) the interference in 
the electron orbital motion via spin-orbit coupling. The crossover between the latter two symmetry classes may be 
achieved by varying the Zeeman splitting. Its features will be described in Chapter [Tv|, together with weak localization 
effects. 

As phase coherence is crucial for the occurrence of interference effects, the influence of symmetry breaking on 
conductance measurements is not observable when the corresponding time scale is larger than the decoherence time, 
t b > T ip- As for finite systems, when decoherence of electrons in a dot or a short wire happens faster than the 



escape into bulk reservoirs, t v < t csc , both the enhanced back-scattering effect and UCFs get suppressed |32 33 
as compared to the conductance quantum, e 2 /h. In these situations, no influence of an in-plane magnetic field on 
quantum transport characteristics of the system would be expected - unless the field affects the decoherence time. 
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This is the case when the latter is caused by spin-flip scattering at paramagnetic impurities |34j . In all foreseeable 
regimes [ |§4] , f35[ , degenerate paramagnetic impurities destroy weak localization via spin- flip scattering, the relevant 
time scale being the spin- flip scattering time r s . Furthermore, they suppress conductance fluctuations []36|-f4(| as 
the flipping of impurity spins effectively changes the realization of disorder in the sample and, thus, leads to a self- 
averaging of the conductance. However, a Zeeman splitting of the magnetic impurity levels by a magnetic field causes 
an energy threshold for electron spin-flip processes which slows down the electron spin relaxation. At a high field, the 
electron spin relaxation is possible only due to very few thermally activated impurity spins, so its rate is reduced by 
the exponential factor, t^ 1 {B) cx t~ x x exp[— g^ lmp ' (j,bB /T], where ^ is the Bohr magneton and g( imp *) the ^-factor 
of the impurities. Thus, a sufficiently strong magnetic field restores UCFs p6| , ^8| -^0| as well as weak localization (if 
the orbital effect of the magnetic field is excluded) . The in-plane magnetoresistance of magnetically contaminated 2D 
semiconductors and other features related to the dynamics of paramagnetic impurities are discussed in Chapter |v| 

For all practical purposes, the study of symmetry breaking rates, determining the parallel field effect on the weak 
localization corrections to the conductivity of a 2D electron gas or a semiconductor wire, can be completed using the 
diagrammatic perturbation theory for disordered systems. Although one can find more details about this technique in 
several textbooks and review articles (see, e.g., J4l|]), for the sake of completeness, we describe its necessary elements in 
Chapter 0. As an alternative, field theoretic methods may be used. In chapter VI, we illustrate the supersymmetric 
non-linear sigma-model in application to one of the problems listed above, i.e., the orbital effect of an in-plane magnetic 
field (via subband mixing) on the transport properties of 2D electrons. 



II. THE DIAGRAMMATIC TECHNIQUE: DIFFUSONS AND COOPERONS 

Diagrammatics is a perturbative method which allows one to a) classify different contributions to the perturbation 
series and b) sum up the relevant terms. To be specific, here - for simplicity - we choose the random disorder potential 
V to be drawn from a Gaussian white noise distribution: 

{V(t))=0, ( V (r)V(r')) = -±-5(r-r'), (9) 

iTtVT 

where v is the density of states (DoS) at the Fermi level and r is the mean free scattering time. Note that throughout 
this review we will use units where H = c = e = I. 

The starting point of the perturbative approach is the representation of the Green function as a series in powers 
of V, i.e., G = Go + Go^^Go) 11 - Under averaging, diagrams with single impurity lines vanish (due to (V) = 0). 
Thus, the diagrammatic expansion involves only diagrams with paired impurity lines. The average Green function is 
then given by a Dyson equation, depicted schematically in Fig. |J: 

oo „ 

(G) =Go + GoV(SG r = Go + GoS(G) <=> (G) = - " 

where the self-energy £ contains all irreducible diagrams, i.e., diagrams that cannot be split into two by cutting one 
Go-line. 
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FIG. 4. Diagrammatic representation of the Dyson equation. 



The dominant contribution to the self-energy reads 

E±(p) = (J dq y(q)G± (p - q; e)U(-q)^ , (10) 

where the shorthand notation dq = d d q/(2Tr) d has been introduced. Since the potential is 5-correlated in real space, 
its Fourier transform does not depend on momentum. Hence, = J dpG^p; e)/(27wr). 

The real part of S, leading to a shift in energy, can be absorbed in the ground state energy. Using that the 
unperturbed Green function reads G^p; e) = (e ± iO + £p) _1 , where £ p = £f — p 2 /(2m), and the definition of the 
density of states, v(e) — L~ d J2 P <5( e + £ P ), the imaginary part of E obtains 
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(11) 



The associated length scale I is the decay length of the average Green function as we will see shortly. In the case of 
weak disorder, t _1 <C ep, all other contributions are small in l/(kp£), and Eq. (^Tj) determines the self-energy £ in 
the self-consistent Born approximation (SCBA). 

Inserting the above result into the Dyson equation yields 



(G ± )(p;e) = 



1 



(12) 



In real space representation, this leads to a decay of the average Green function on the scale of the mean free path, 

(G)(r,r') = G (r,r')e-^. (13) 

Having found the average Green function, the next step is to calculate correlation functions of the form[j] 

F(pi,p , 1 , P 2,p , 2 ;a;) = (tT(pi, p[;e)G+ (p 2 ,p' 2 ;e + u)). (14) 

Again - as for the averaged Green function - the diagrammatic perturbation series involves summing up diagrams 
with paired impurity lines. The dominant contributions are series of ladder diagrams and maximally crossed diagrams, 
see Fig. g. With the notation (AB} C — (AB) — (A)(B), i.e., subtracting the disconnected part of the correlator, the 
contribution of connected diagrams can be written as 

(G-(p 1 ,p' 1 ;e)G + (p 2 ,p^; £ + ^)) c (15) 
= (G-)( P1 , e) (G-)(pi, e) (G+)(p 2 , e + w) (G+)(p 2 , e + «))r(p 1>Pai pi, p 2 ; u) 5(pi - p 2 - pi + p 2 ), 

thus defining the reducible^ vertex function T(pi, p 2 , pi, p 2 ; to). 



p + q 


p'+q 


p + q p"+ q 


p'+q 


p + q p"+q 


p"'+q 


p'+q 








+ 




+ 


+ . . . 


P 


P' 


P P" 


P' 


P P" 


P'" 


[>' 




-P'+Q 


-P + Q 


-p'+Q -p"+Q 


-p + Q 


-P'+Q -p"'+Q 


-p"+ Q 


-P + Q 








+ V 




+ ^ 






+ . . . 


















P 


< 

p' 


p p" 


p' 


P p" 


P'" 


p' 





FIG. 5. Ladder and maximally crossed diagrams. 

Consider first the sum of ladder diagrams depicted in the upper part of Fig. Due to momentum conservation 
at each vertex, the momentum difference q = Pi— p 2 (= pi~ p 2 ) is constant, and T depends on this difference only. 
Then T(q, lu) is given by a Bethe-Salpeter equation, the two-particle analogue of the Dyson equation: 

r(q, w) = 7 + 7 [ / dp" (G~)(p" + q, e + w) (G+) (p", e + w) 1 T(q, «), 



n(q,«) 



The connected part of averages of the form (G + G + ) or (G~G~) vanishes, i.e., no long-ranged correlations exist. 
2 In the case of two-particle functions, a diagram is called reducible if it can be split into two separate diagrams by cutting 
two (G)-lines. 
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where 7 = \/(2itvt). Anticipating the result that T diverges for q, u> — > 0, we can approximate the irreducible vertex 
function II for small |q|,u>: 

IT(q, u>) ~ 2ttvt(1 + iujt - Dq 2 r), 

where D = v\t j d is the diffusion constant (d dimensionality of the system). Thus, the sum of ladder diagrams finally 
yields 

^ = 1 777 7 = o 1 2 n 2 1 ■ = w )- ( 16 ) 

1 — 7ll(q, lj) 2-kvt 1 Dq* — iuj 

This is a diffusion pole or so-called diffuson. It can easily be seen that the diffuson is not affected by a weak magnetic 
field: Minimal coupling implies that the magnetic field shifts all momenta by the corresponding vector potential, 
p -> p - A; however, (pi - A) - (p 2 - A) = pi - p 2 = q. 

A maximally crossed diagram Q (for a detailed discussion see e.g. Q) can be converted to a ladder by reversing 
one G-line. Since now all the arrows point in the same direction, momentum conservation at the vertices requires the 
momentum sum Q = P1+P2 to be constant. Apart from that, the structure of all equations is the same as for the 
ladder diagrams. Thus, 

This is called a Cooperon - in analogy to superconductivity, as it corresponds to a correlator in the particle-particle 
channel (as opposed to the particle- hole channel for the diffuson). Now the presence of a magnetic field requires 
replacing Q — > Q — 2A which leads to a decaying of the Cooperon. 

Taking into account the spin of the electrons, diffusons and Cooperons split into singlet and triplet modes. This will 
become important when considering spin-orbit coupling or spin scattering. Details will be given in the corresponding 
chapters |[y| and |y|. 

In coordinate representation, the expression obtained in Eq. (|l7j ) has the form of a diffusion equation, 

[£>(V + i2A) 2 -iuj]C = (S(x-x'). (18) 

(A similar equation holds for the diffuson, Eq. (|l6|).) In a finite size geometry, this equation should be complemented 
by boundary conditions. For an 'open' surface, as, e.g., the interface with a bulk reservoir, the boundary condition 
reads 

C = 0. (19) 

For an insulating surface, for example, the edge of a metallic dot, the corresponding boundary condition expresses 
the lack of a phase-coherent current density across the sample surface, 

ny • (V + i2A)C = 0, (20) 

where ri|| is a unit vector normal to the sample surface. In the case of zero magnetic field, this boundary condition 
indicates that the lowest Cooperon and diffuson modes correspond to the constant (Q — 0) solution and are gapless. 
This choice of the lowest mode, being constant in space, is called the zero-dimensional (0D) approximation. In the 
presence of a magnetic field, one has to give the analysis of the lowest modes in closed systems more attention, applying 
a gauge transformation to the gauge where ny - A = 0, and only then making the 0D approximation. Subsequently, 
the discrete spectrum of relaxational modes describing the Cooperon and diffuson can be found using the standard 
Hamiltonian perturbation theory methods. This is a useful trick which enables us find the Cooperon and diffuson 
functions and, in particular, their lowest eigenmodes in finite size geometries. 



III. IN-PLANE MAGNETORESISTANCE EFFECT DUE TO INTERSUBBAND MIXING 

As pointed out in the introduction, the influence of an in-plane magnetic field on the quantum transport charac- 
teristics of a 2D or quasi-2D system can be divided into orbital and spin-related effects. In this section, we consider 
the orbital effect - treating electrons as spinlcss. 
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There are two features that make the issue of the influence of an in-plane magnetic field on quantum interference 
effects in inversion layers and thin metallic films non-trivial. The first issue is related to the fact that a finite 
spatial extent of electron wave functions along the confinement direction, i.e., the z-direction perpendicular to the 
2D plane, is needed in order that the electron can accumulate an additional magnetic field-induced phase if the field 
is exactly parallel to the 2D plane. Furthermore, the possible accumulation of such an Aharonov-Bohm phase is 
sensitive to the inversion symmetry properties of the system. In particular, for a quantum well with a potential 
profile that is symmetric with respect to the transformation V z : z — ► —z, an in-plane magnetic field cannot induce 
any dephazing, and the field-dependent rate r^ 1 is equal to zero. This property reflects the Berry-Robnik symmetry 
effect [|4|: Consider a system described by a Hamiltonian H that at B = is invariant under a discrete symmetry 
transformation, such as V z , and at finite fields remains invariant under the combined time and z-coordinate inversion 
transformation, V Z T , i.e., V Z T'H(B) — H(B). Then there always exist pairs of paths having exactly the same value of 
the semiclassical action which, therefore, interfere constructively even in the presence of a magnetic field. As a result, 
the suppression of weak localization by an in-plane magnetic field is possible only if either the scattering potential in 
the quantum well is z-dependent |l3| ] (since a generic z-dependence does not respect Pz-symmetry) , or if the form of 
the confinement potential in z-direction has no inversion symmetry Jl5|-p^7t (as usually is the case for a inversion layer 
in a heterostructure). In the first subsection, we provide a quantitative analysis of this effect. In the second subsection, 
this analysis is extended to multi-subband systems Ji5| , ^6[ , and the results are compared to the magnetoresistance in 
thin films with diffusive transverse motion. 

The second non-trivial issue concerns the so-called geometrical flux cancellation characteristic for ballistic metallic 
films (wires), where electrons scatter only from surface (edge) defects at two parallel surfaces. In such a 'quasi-ballistic' 
film, the finite resistance is provided by the diffusive scattering of carriers from the surface. Quantum corrections 
to the conductivity are related to the interference between waves propagating along paths consisting of a sequence 
of ballistic flights between the film edges. Among those paths, the trajectories containing ballistic segments that 
are almost parallel to the film surface and, therefore, much longer than the film width d play an important role. 
These so-called Levy flights (in the terminology of anomalous diffusion theory jf5|) produce a logarithmically large 
contribution to the sample conductance [ff6|. The oriented area encircled by a closed loop made of quasi-ballistic 
paths is exactly equal to zero thus suppressing the magnetic field effect on the interference of electron waves 
propagating along them. An electron can accumulate a magnetic flux only due to a curving of its trajectory by the 
magnetic field itself. This makes the role of Levy flights important for the interference-related magnetoresistance of 
ballistic films and wires: curving of the longest ballistic flights - which is a purely classical effect - switches their 
scattering at the surface roughness from one film surface (or wire edge) to the other p5j . This peculiar situation 
and the resulting unexpected influence of a magnetic field on the localization properties in quasi-ballistic wires are 
discussed at the end of this section. 



A. Phase-breaking by an in-plane magnetic field in 2D inversion layers 

The influence of an in-plane magnetic field on the orbital motion of carriers in a heterostructure or quantum well 
is a result of the finite width d of any 2D layer. This effect has been discussed previously in various contexts [p~8|— ^T]| . 
The Lorentz force on the electrons generated by the in-plane field Bp mixes up the electron motion along and across 
the confinement direction, thus resulting in a modification of the 2D dispersion, E(p) — > S(Bii,p). In particular, 
a) the 2D electron mass increases in the direction perpendicular to B||. Furthermore, in heterostructures with a 
confining potential that does not possess inversion symmetry, b) the in-plane field also lifts the p — > — p symmetry in 
the dispersion law M,tt9[: 

£(B||,p)-E(B||,-p)oc (p-[B|| xe 2 ])Vo, 

where e z is a unit vector perpendicular to the plane. This leads to a field-dependent suppression of the WL signal as will 
be discussed below. By contrast, for the z-inversion symmetric problem, due to the Berry-Robnik phenomenon [E|, 
the Cooperon mode remains massless, and there is no magnetoresistance. 

In the presence of an in-plane magnetic field, the effective 2D Hamiltonian for electrons in a heterostructure with 
the confining potential profile W{z) can be obtained from the 3D Hamiltonian, 

— A) 2 f) 2 

H 3D = ^^-^ + W(z) + u(r,z), (21) 

using a plane wave representation, = e tpr (f>^ (z) for the electrons in the lowest subband. Here, 4>^\z) are the 
eigenfunctions of the z-dependent part of the Hamiltonian. In particular, at B\\ = 0, the in-plane and perpendicular 
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motion are separable, and the corresponding eigenfunctions are given by |fc) ee 4 0) . Due to the mixing between 
subbands |0) and \k >0) by an in-plane magnetic field, the eigenfunction (f>^ (z) in the presence of the field depend on 
the in-plane momentum p. In order to find (^^{z) and the corresponding energy i?(By,p) for each plane wave state 
a perturbation theory analysis - as described in the following - or a numerical self-consistent-field technique may be 
used. 

Furthermore, in Eq. (|2l|), the vector potential reads A = (z — 2o)By x e z , where zq = (0|z|0) is the center of mass 
position of the electron wave function in the lowest subband (at By = 0), and the potential it(r, z) is a combination 
of the Coulomb potential of impurities and the lateral potential forming the quantum dot. 

For weak to intermediate magnetic fields By, the effective 2D Hamiltonian takes the form 

^ 2D = hi ~ P ^ 7(B|1 } + } + u(r) - (22) 

Here, p = — iV — a(r) is a purely 2D momentum operator and p± = p • x e z ] its component perpendicular to By. 

Thus, the application of an in-plane magnetic fields results in a) an effective vector potential a(r), and b) two 
additional terms in the free electron dispersion due to p^-dependent inter-subband mixing. The first term, STt^ — 
— p5_7, lifts the rotational symmetry by causing an anisotropic mass enhancement Jl8| , ^9[ . It increases the 2D density 
of states and, for a 2D gas with fixed sheet density, it reduces the Fermi energy calculated from the bottom of the 
2D conduction band, i.e., Ep(B») = Ep — j(Bii)p^/2. The second term, SH^ = Pj_f3, is related to the time-reversal 
symmetry breaking by By . A perturbative analysis of the problem results in the field dependences 7 ~ B?, and 
(3 B^. Indeed, an expansion of the plane- wave energy up to the third order in p± yields 



k> 1 



Rin ^ S ll V (0\z\k)(k\z\k')(k'\z\0) d ( d 



m3 k W>\ ( £ fc ~ £ o)(efe' - £o) m \Xi 

where are the subband energies and Xb — 1/ \/B\\ is the magnetic length. 

In TI2D given by Eq. (|22j), the "in-plane" disorder is incorporated in the form of a scattering potential u(r) w 
(0|u(r, z)|0). It can be characterized by the mean free path, i 3> 1/pf or, equivalently, the momentum relaxation time 
r, related to the diffusion coefficient D = Vpr/2. The modification of the electron density of states by By only slightly 
affects the value of the electron mean free path. However, the presence of a parallel field also changes the symmetry 
of the Born scattering amplitudes between plane waves, / pp ' = (\l/p|u(r, z)|\& P '). Due to the momentum-dependent 



subband mixing, / pp < acquires an additional contribution, / pp ' = / pp , {l + {pi_ + p' ± )B\\(} , where 



c = L - y 



_ AO) 
- /pp' 

(0\u(p- P ',z)\k)(k\z\Q) 



u(p-p',z)\0) ^ tfc-eo 



This corresponds to the presence of a random gauge field in the effective 2D Hamiltonian |l3],|l4|,[l7| , 

(0\u(r,z)\k)(k\z\0) m p 

which can be interpreted as a result of an effective 'curving' of the 2D plane by impurities with a z-dependent 
scattering potential in the system. In the presence of an in-plane magnetic field this curving generates a random 
effective perpendicular field component, b± = [V x a] . In systems, where scattering is dominated by Coulomb 
centers separated from the 2D plane by a spacer and, thus, is almost independent of z, a smaller effect may be taken 
into account, namely Sa ~ [By x e z ]([By x e z ] ■ V) 2 u(r). However, 5a has a negligible influence on the quantum 
transport characteristics of 2D electrons as compared to the effect of the field-dependent electron dispersion. 

The perturbative calculation of two-particle correlation functions, i.e., diffusons 2?(x, x!;u>) and Cooperons 
C(x, x';cj) requires taking into account all field-dependent terms in the effective 2D Hamiltonian Eq. (0). 

Admitting for different magnetic fields By! and Bp, the result for the diffuson has the form 

[-UV 2 -iu + T^]V = <5(x-x'), (23) 
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where Co = to + 8 with 5 = pp[y(Bi\i) — 7(-Bi|2)]/2. The field-dependent rate T d 1 , B\\2) reads 

*d _1 = ^tiPp) ~ 7(S|| 2 )] 2 + £p|[-B yi - B|| 2 ] 2 (24) 

which vanishes for B^ = Bn 2 . 

The first term in Eq. ( |24] ) comes from a deformation of the Fermi circle by the magnetic field B\\ . The second term 
takes into account the field effect on the scattering of plane waves. Furthermore, Eq. ( p3| ) contains the difference 
between the electron kinetic energies in two measurements of conductance, Co = Ep(B^i) — _Ef(-B|| 2 ), each of them 
shifted by the magnetic field with respect to the Fermi energy Ep at B\\ = of the electron gas with the same sheet 
density. The latter fact is important, since, for lateral dots where the electron density is fixed, one should substitute 
uj = 0. 

The Cooperon equation can be represented in the form 

[D(-iV - A ) 2 + 1 + r" 1 ] C = <5(x - x'). (25) 
It contains an additional decay rate t~ 1 (B^i, -Bp) given as 



1 _ t Pf \ f3(B u ) + P(B n2 ) 
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C z p z f 

2r 



[B\\i+B\\2] ( 26 ) 



which accounts for the dcphazing of electrons encircling the same chaotic trajectory in opposite directions. Thus, the 
result of lifting time-reversal symmetry by an in-plane magnetic field can be described by the Cooperon phase-breaking 
rate 

r-,f = T-\B h B«) = ^f3 2 (Bn) + 2 -^B\ ~ aB\ + bB\. (27) 



Note the gauge shift in the Cooperon equation (|2^). The 'vector potential' A = |p|m/3(B||)[-g|j- x e 2 ] results from 
the following artifact: the cubic term in the effective electron dispersion not only lifts the p — > — p inversion symmetry 
of the line E(B^ , p) = Ep(B^), but also shifts its geometrical center with respect to the bottom of the 2D conduction 
band. However, since for the conductance only electrons with energy e = E-p matter, such a shift can be eliminated 
by choosing a slightly modified initial gauge. This can be corrected for now by applying a gauge transformation 
C — * e lqx C directly to the Cooperon. In other words, the phase-coherent transport is only affected by the B\\ -induced 
p — > — p asymmetric distortion of the Fermi circle into an oval, but not by a shift of its center in the momentum space. 
That is, only the -By -dependent anisotropy of the electron wavelength along the Fermi line affects the interference 
pattern of current carriers. Note that both terms in Eq. ( p7j ) are non- vanishing only in the absence of z — > — z inversion 
symmetry^] 

The phase-breaking rate in Eq. (p7| ) can also be used to describe magnetoresistance and universal conductance 
fluctuations (UCF) in lateral semiconductor dots. For quantum dots, weak localization corrections and the variance 
of UCF can be represented as 



Swl(B||) cx J (dx) W(x)C(x,x;0) = (ff(B|,)) - (g) u , 
where the subscript 'it' stands for unitary, and 



(Sg(B n )Sg(B ]]2 )) cx / \dx)(dx') W "(x)W (x') ^ |^(x,x';c)| 

•J t~> t-\ r> 



v=v,c 



where lo — Ef(B»i) — Ef(B\\ 2 ), and Ef(B\\) is the Fermi energy of the 2D electron gas calculated from the bottom 
of the conduction band. The dispersionless weight factors W(x) take care of the particle number conservation upon 
diffusion inside the dot Ej| and incorporate the coupling parameters to the leads. In the zero-dimensional limit, both, 



3 In fact, the first term requires an asymmetric confining potential whereas for the second term a generic z-dependence of the 
impurity potential is sufficient. 
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gwL and {6g(Bi\i)Sg(Bp)} , are dominated by the lowest (spatially homogeneous) Cooperon (diffuson) relaxation mode 
Ao, which is determined by the rate of escape to the reservoirs, r~* <C D/L 2 (Thouless energy), and the Cooperon 
suppression by time-reversal symmetry breaking described by t^m . The latter expresses the efficient reduction of the 
fundamental symmetry of the system from orthogonal (o) to unitary (u). In the absence of spin-orbit effects, the 
parameters r~* and describe the value of WL corrections as well as the variance of UCF, (Sg 2 (B^)), as compared 

to their nominal values, <7wl(0) = (g) ~ (g)u and (5g 2 ) u : 



9wl(B\\) = ,9wl(0) — — , 

1 + t csc /t b 

(6g 2 (B {] )) = {Sg 2 }Jl + —± 



(28) 
(29) 



Thus, these two parameters can be studied from the UCF fingerprints measured by changing the shape of a dot in 
multi-gate devices or by varying the Fermi energy in back-gated dots. 




1 1.5 2 2.5 

n(10"cm" 2 ) 

FIG. 6. Calculated dependence of the parameters 7 and j3 on the sheet density of 2D electrons. The insets show the effect 
of Bu on the symmetric, [E(p) + E(—p)]/2, and anti-symmetric, [E(p) — E(—p)]/2, parts of the 2D electron dispersion in a 
broader range of Bu , where a perturbative expansion is not applicable. 



In ballistic billiards or heterostructures with a z-independent impurity scattering potential, the in-plane magnetic 
field effects are governed by the unusual B^-dependence of t^I p5^p7t that we attribute to the effect of a cubic term 
in the 2D electron dispersion generated by the field. This should give rise to a relatively sharp crossover between 
'flat' regions related to regimes with orthogonal and unitary symmetry, respectively. To get an idea about the 
relevance of this effect for the large-area (~8/xm 2 ) quantum dots studied by Folk et al ||, quantitative estimates for 
the parameters 7(-Bji) and /3(-B||) can be obtained by evaluating the electron dispersion within a fully self-consistent 
numerical method |L7|. Using the nominal growth parameters of the Al.34Ga.66As/GaAs sample studied in Q, the 
dependences of 7 and on the in-plane magnetic field for an electron sheet density of 2xlO n cm~ 2 shown in the 
insets of Fig. || a) and b) were found. At low fields, 7 ~ B 2 and f3 ~ £?jj , as anticipated in the perturbation theory 
treatment. The proportionality coefficients are plotted in Fig. ^ a) and b) versus the electron sheet density. Both, 
the effective mass renormalization in the quadratic term of the energy dispersion and the time-reversal symmetry 
breaking cubic term are larger at lower 2D electron gas densities, due to the weaker confining electric field, i.e., the 
increased width d of the potential well. From this analysis, one estimates the field necessary to suppress the weak 
localization effect completely as Bu = 0.6 ~ 0.8T, which is in agreement with the observed complete suppression of 
WL corrections beyond B\\ ~ 0.6T in experiment fl49| . 
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B. In-plane magnetoresistance in inversion layers with several filled subbands and in diffusive metallic films. 

If the width of the quantum well exceeds the Fermi wavelength of the electrons, several subbands become occupied 
- as shown in Fig. ^ - and, thus, contribute to the lateral transport. As for the one-subband case discussed above, 
the effect of an in-plane magnetic field depends sensitively on the microscopic structure of the wavefunctions in the 
z-direction perpendicular to the plane, that determines possible couplings between the subbands. In particular, in 
the absence of z-dependent impurity scattering, the influence of By is different for quantum wells with a symmetric, 
W{z) = W(—z), and asymmetric, W(z) 7^ W(—z), confinement potential due to the Berry- Robnik symmetry effect. 




FIG. 7. Schematic picture of the quantum well. Two exemplary subband wavefunctions are shown. The spatial profile of 
the impurity potential is sketched on the bottom of the well. 



In the absence of an external magnetic field (Bt\ = 0), the subbands are uncoupled and, thus, each subband 
contributes separately to the conductivity. The weak localization correction for a film with M occupied subbands is, 
therefore, given as 

e 2 

g W h = M — \ti(t / Tfi) , 
irn 

where phase coherent propagation in each subband is limited by the inelastic decoherence time r^ 1 oc T p . 

An external magnetic field now plays two complementary roles: a) it breaks time-reversal (T) symmetry and b) it 
couples the different subbands (i.e., if the random disorder potential is z-independent, it represents the only possible 
coupling between subbands). In fact, these two roles are linked: if the subbands remained completely decoupled, 
T-invariance of the 2D motion in each of the M subbands would be preserved because the vector potential of the 
parallel field can be gauged out in each particular subband. The field-induced subband mixing results in a spectrum of 
effective dephazing rates 1/tq (where fc = 0, . . . ,M— 1) with the following properties: at least M— 1 dephazing rates 
are quadratic in the magnetic field, \/t^ >0 ^ oc B 2 . The remaining is exactly equal to zero for a 'Pz-symmetric 

system whereas it acquires a quadratic field dependence as well if the symmetry is broken (due to an asymmetric 
confining potential or z-dependent impurity scattering). As a result, the field and temperature dependence of the 
magnetoconductivity are different for a symmetric vs asymmetric confinement potential, 

e 2 ( plnT + 2(M-l)lnB l{ for W(z) = W(-z), 
gwL (B h T)- — ^ 2MlnB|| for W(z) ? W(-z). (30) 

The derivation of the above weak localization magnetoconductivity formula is most conveniently performed within 



the non- linear sigma-model formulation that will be presented in chapter VI. 

For strong inter-subband impurity scattering, i.e., when the inter-subband scattering rate exceeds the subband 
splitting at the Fermi level, the weak localization suppression by the in-plane magnetic field happens in the same 
way as in thin metallic films with diffusive transverse motion, d ^> I. In a film with diffusive transverse motion, the 
electron encircles a random flux of order $ ~ B\\d 2 during its diffusive transverse flight, where the time of flight is 
given as t z ~ d 2 /D. While diffusing along the wire for a longer time, t 3> t z , the electron encircles a random flux 
with the variance 

($2) „± x (B ]] d 2 ) 2 =tD(B l] d) 2 , 

7~z 
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which results in a field-induced phase breaking at the rate |22| 

Ts 1 - D(B {] df. (31) 

Note that the same geometrical argument is applicable to a thin diffusive wire with cross sectional dimensions ~ d in 
a magnetic field applied along the wire axis 




FIG. 8. Relevant areas for perpendicular and parallel magnetic fields. This admits for estimating field dependence of the 
respective phase breaking rates, 1/tb ± ~ l/'x ~ B± and 1/ts,, ~ ~ By. 



Quantitatively, tb can be found from an analysis of the lowest eigenvalue of the Cooperon propagator. The Cooperon 
propagator is defined by the (diffusion) equation 



with the boundary condition 



D{iV r - 2A(r)) 2 + — \ C(r, r') = S(r - r') 

T <t>. 



(V n + 2iA-n)C = 0. 



Choosing the gauge A = Buze y , one can separate variables in the eigenfunction equation for the Cooperon modes, 



-D 



o2 

Q^~Ql-(Qy-2B llZ ) 2 



(f>k,Q(z) = ^fe,Q0fe,Q(^), 



and find the lowest eigenvalue f^o.Q = DQ 2 + t b 1 within perturbation theory. This yields the magnetic decoherence 
rate H 

= ±D(B l{ df 

which can be substituted into the weak localization formula for a diffusive film. A similar calculation gives t^ 1 — 
D(B\\d) 2 /8 for a circular diffusive wire with diameter d. 

C. In-plane magnetoresistance in pure metallic films and quantum magnetoresistance in quasi-ballistic wires 

with rough edges. 



For a clean metallic film or wire, where the dominant scattering processes involve surface or edge roughness, the 
estimate of tb given in Eq. ( pl| ) is incorrect due to the exact geometrical phase cancellation described in p3[ . Thus, a 
magnetic field Bn cannot affect the interference pattern, i.e., neither WL magnetoresistance nor magneto- fingerprints 
are to be expected in these systems at any orientation of the magnetic field. As a result, magnetic field effects acquire 
several unusual features. 

Here, we consider a quasi-ballistic wire with width d and length L. In such a wire, a diffusive path may contain 
ballistic segments with length rj = d/ sin a much longer than the typical length scale d, crossing the wire at the small 
angle a. These segments appear with the probability 



dr/(a) 



da 



V 



for L > 1] 3> d Xp. Thus, the mean free path t in such a wire exceeds the width of the wire by a large logarithmic 
factor, £~ (j]} ~d \n(L/d). 

Quantum in-plane magnetoresistance in such a wire may appear only due to a curving of the longest ballistic 
paths by the field itself. The effect of this curving leads to a switching of the longest segments between the edges 
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of the wire (or surfaces of a film) and, thus, to a sudden suppression of the interference contribution from paths 
including these longest ballistic flights. However, the change in the weak localization part of the conductivity in such 
a situation is overcome by the change in the classical part of the conductivity: due to the cyclotron curving of free 
electron trajectories in combination with the finite width of the wire, the longest ballistic flights are suppressed, i.e., 
i] < yjr c d - yielding a mean free path t{B\\) ~ hdln(r c /d), where r c = pp/Bn. Thus, the shortening of the longest 
ballistic flights happens simultaneously with their involvement into phase accumulation, immediately bringing in a 
large time-irreversible phase of order of the total flux penetrating through the entire sample. On the other hand, 
paths that encircle a large flux can participate in the formation of a random interference pattern. Therefore, starting 
from the field value -Bdefl = <f>od / '{Li 2 : Af) (implying y/r c d < L) - that coincides with the set in of a classical (positive) 
magnetoresistance - the quantum conductance should acquire a universal random dependence on the field with a 
characteristic scale related to a flux change of order of the flux quantum. By the same reason, the magnetic field 
effect on the localization properties of quasi-ballistic wires is opposite to the commonly expected crossover between 
orthogonal and unitary symmetry classes, i.e., the magnetic field tends to shorten the localization length in a quasi- 
ballistic wire. 

These qualitative expectations have been verified numerically in Ref . . The results of Ref . pH] are based upon 
the numerical solution of a two-dimensional Anderson Hamiltonian on a square lattice, H = J^. l*) £ i(*l l*)(j |> 

where (ij) denotes nearest neighbor sites i and j. The structure considered consists of two ideal leads attached to a 
scattering region that is W sites wide and L sites long (all lengths are in units of the lattice constant a). For simulating 
bulk disorder, the energies 6j in the scattering region are taken uniformly from the interval — U/2 < &j — eo < U/2, 
where U is the disorder strength. For sites on the boundary, e, = eo + £ s with es — 10 4 . The rough structure of the 
boundary was generated by having an equal probability of either 0, 1 or 2 sites at each edge with on-site potential 
£q + £b- To simulate the effect of a magnetic field, a Peierls phase factor has been incorporated into V in the scattering 
region. 
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FIG. 9. (a) Distribution function P(L/£) for 4 quasi-ballistic samples with W =15 and L =52(A), 104(B), 208(C) and 416(D) 
and for a wire (U) with on-site bulk disorder U — 2.0 (nominal mean free path £ « 8.5, W = 15 and L = 52). (b) Distribution 
function P(g) for the two quasi-ballistic structures A and C. (c) E^{L) for quasi-ballistic wires with eo G [1.5, 1.7] (circles) 
and eo £ [1-0, 1.2] (diamonds) and for a wire from the U series (squares). Also shown is the analytical result (cf. Eq. ( |32j ) with 
hv? = 12.3 and W — 9.3 (solid curve). The dashed curve corresponds the ergodic law E Th ~ hD/L 2 with hD = 62.5. 



The effect of Levy flights has been identified in the numerical results in several ways. One was to study the 
distribution of transmission coefficients through the wire and the distribution -P(£ _1 ) of the corresponding Lyapunov 
exponents This is shown in Fig. ^(a) (taken from Ref. p^ ]) for 4 series of quasi-ballistic wires (A-D) as well as a 
series of samples with bulk 'defects' (U). As pointed out by Tesanovich et al |5(| and verified numerically in Ref. fl5l| , 
the length of Levy flights in quantum systems is limited. The limitation is due to the fact that the uncertainty in 
the transverse momentum, 5k± ~ in a wire with a finite width sets a quantum limit to the angles a ~ d/rj 
which can be assigned to a classically defined ballistic segment, namely a > 8a ~ Sk±/kp ~ Xp/d. This sets the 
cut-off ?7 max = d 2 j 'Ap and entails a finite localization length L c ~ 77 max - Samples from the series A and B meet the 
criterion L < L c which manifests itself in the distribution of Lyapunov exponents by a pronounced peak at small 
£ -1 corresponding to eigenvalues T n i=s 1. This should be compared to the plateau-like distribution [j52| obtained in a 
sample from the series U using the same numerical procedure. Note the finite width of the ballistic peak in -P(£ _1 ) 
and that P(0) = 0. The enhanced density of small can be identified even in samples from the series C and D 
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with L w L c . In these samples, -P(£ _1 ) starts to show a periodic modulation specific to the localized regime, where 
the spectrum of Lyapunov exponents tends to crystallize [5^j5^ | . 

The distribution of eigenvalues of the transmission matrix results in a finite- width distribution of conductances and, 
hence, conductance fluctuations. The statistics of conductance fluctuations, P(g), for wires with edge disorder is shown 
in Fig. ||(b): the results of numerical simulations p5| are plotted for series A, where (g) » 2.0 and the distribution 
is almost Gaussian, and series C, where (g) » 0.71. (The conductance is expressed in quantum units e 2 /h.) The 
distribution function P(g) is the result of the analysis of various realizations of disorder. When calculating correlation 
functions, there is an additional averaging over energy, since the energy dependence of conductance fluctuations 
is random on the scale of the Thoulcss energy Et\i- The latter can be determined from the wire conductance as 
g ~ ETh(L)/A(L), where A is the mean level spacing in the wire. Taking into account the logarithmic multiplier in 
the dependence of the classical conductance on sample length due to the Levy flights, one obtains 

This result for the wire conductance implies an anomalous scaling of the Thouless energy, i?Th = irdvpL - 2 ln(y4b), 
with the length of the wire. In Fig. |9|(c), it is shown that, indeed, the Levy flights manifest themselves in the correlation 
energy of UCF in the quasi-ballistic regime. For comparison, it is demonstrated that the results for bulk-disordered 
samples (U) coincide with the standard scaling law Et^L) ~ L~ 2 . 
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FIG. 10. (a) Conductance for a quasi-ballistic wire (solid curve) from the B series and a disordered wire (dashed curve) with 
mean free path I ~ 16. (b) Variance var(<?(_B)) for quasi-ballistic wires from the B series and disordered wires. 



Coming now to the magnetic field effects. The expected deflection of Levy flights at B > i?dcfl entailing a break 
of geometrical flux cancellation is observable in Fig. |l^ as pronounced magnetoconductance fluctuations Sg(B) = 
g(B) — (g) beginning together with a gradual decrease of the average conductance value - which is a purely classical 
effect. Their variance corresponds to the usual UCF value, and the auto-correlation function as a function of magnetic 
field gives the correlation field B c corresponding to a change of magnetic flux through the sample area of order 2.5</>o 
as compared to 1.5</>o in the bulk-disordered case. This seems to explain an earlier experimental observation of UCF 
in quasi-ballistic semiconductor wires f35f . 




B (t/a 2 ) V% B (^a 2 ) 

FIG. 11. (a) The ensemble averaged conductance as a function of magnetic field for wires from the A and B series, (b) The 
same curves (A-D) as in Fig. [](b), but at finite magnetic field B = 0.02</>o/ 'a 2 . (c) The localization length for quasi-ballistic 
wires (solid curve) and disordered wires (dashed curve) with mean free path t » 8.5. 



Furthermore, Fig. [Ll] shows the tendency of the localization length in a quasi-ballistic wire to increase with the 
magnetic field - as opposed to the behavior of L C (B) in systems with diffusive transverse motion |4l]j5^| . The relevant 
field scale at which the localization length in quasi-ballistic wires starts to shorten, B c ~ 0oAf/IF 3 , can be obtained 
from -Bdcfl by replacing the sample length L with L c . Formally, this yields a similar field scale as the one leading to 
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a crossover in L c between two symmetry classes due Aharonov-Bohm type phase effects. Thus, the latter is hindered 
by the deflection effect and does not appear as an intermediate crossover regime. 



IV. IN-PLANE MAGNETORESISTANCE DUE TO SPIN-ORBIT COUPLING 

A. Spin-orbit coupling in 2D heterostructures and quantum dots 

Coupling between the electron spin and its orbital motion is a relativistic effect inherent to many metallic and 
semiconductor systems. For example, in zinc-blend type III-V semiconductors, such as GaAs, which have no inversion 
symmetry in the unit cell, the 3D bulk dispersion of conduction band electrons contains terms that are cubic in 
the electron momentum and linear in its spin operator ]5(| . This generates a linear spin-orbit coupling for electrons 
confined to the 2D plane in a heterostructure or a quantum well 57 Bl| , 



H ° = 2^ + 2^ [PX ^-v + ^iP^-Py-y) = ^ + •£> (33) 

where n p characterizes the direction and frequency of spin precession (of an electron moving with momentum p) 
measured in units of uf/Aso, where 1/Aso = y "Ir + o^.. 

The form of the spin-orbit (SO) coupling in Eq. (Q) is specified for a GaAs heterostructure grown along the 
(001) plane which has the symmetry of a square lattice without inversion symmetry within the unit cell. It includes 
two possible combinations of spin and momentum operators that are invariant under the symmetry transformation 
corresponding to such a lattice symmetry, parameterized by two constants, a B R (the so-called Bychkov-Rashba term) 
and a CT (the crystalline anisotropy term) . 

For a particle moving diffusively in an infinite heterostructure, spin-orbit coupling provides an efficient mechanism 
for spin relaxation. In the frame of the moving electron, its spin undergoes a precession with angular frequency 
£2 ~ npfp/Aso which randomly changes direction each time the particle changes its direction of propagation due to 
scattering at impurities. Using the parameterization of the SO coupling by a spin-orbit length Aso> as in Eq. (|33|), 
the precession angle acquired between successive scattering events can be estimated as 5<p so ~ tvf/\so = ^/Aso- 
For weak SO coupling, £/\so "C 1, the electron spin precession can be viewed as a random walk of the electron spin 
polarization vector around the unit sphere with a typical step size <5<^ so . This causes a loss of spin memory at the 
(Dyakonov-Perel) rate [^6| 

^dp ~ D/\l Q . (34) 

In a 2D electron gas, this parameter limits the diffusion time in the triplet Cooperon channel (whereas the singlet 
Cooperon channel is unaffected). Therefore, it suppresses the triplet contribution to the weak localization correction 
formula and controls the crossover between the weak localization (at weak SO coupling) and weak anti-localization 
(at strong SO coupling) regimes in the quantum correction to the conductivity [ p7| , p8[ . 

When the electron motion within the 2D plane is also confined laterally, the use of the so called Dyakonov-Perel 
formula, Eq. ([Mj), dramatically overestimates the rate of spin relaxation. In particular, if the electron motion is 
confined to a single ID channel (i.e., a quantum wire), no real spin relaxation takes place - despite the electron spin 
precession. In this case, the spin precession is a reversible process that rotates the initial electron spin to the same 
definite state at each point along the wire, no matter how many times the electron moves forwards and backwards. 
Mathematically, the above-mentioned property of the SO coupling becomes transparent in the following way: Consider 
a one-dimensional electron moving, say, along the cc-direction. Applying the unitary transformation, 

/ i£ 

U = exp I ^— n(l x ) • a 

where 1^ is a unit vector along the wire axis, the SO coupling is eliminated completely from the electron dispersion 
represented in a rotated spin-coordinate frame J5^,^0|. Therefore, it does not affect any charge transport properties. 
By contrast, the Dyakonov-Perel relaxation of the spin of a 2D electron (with the rate described in Eq. ([54])) is the 
result of the diversity of paths the electron may use to move between two points in the sample. 

A small quantum dot with dimensions L <C Aso represents an intermediate situation between the ID and the 2D 
system as described above. Here, we shall consider a dot coupled to metallic leads via two contacts, left (I) and right 
(r), each with N^ r > 1 open orbital channels. Then the escape rate [p0| from the dot into the leads is given as 



1G 
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(Ni + N r )A/{2ir), 



where A = 2n/(mA) is the mean level spacing (A area of the dot). 

Although there is no unitary transformation that eliminates the SO coupling completely by changing to a locally 
rotated spin coordinate system, the linear term in the small SO coupling constant can be cancelled by a proper spin- 
dependent gauge transformation fl3l[ , thus, generating only contributions to the free electron Hamiltonian of order 
AgQ (or higher). As a result, the influence of spin-orbit coupling on the transport characteristics of a small dot is 
suppressed as compared to the infinite 2D system with similar material parameters. 

However, the application of an in-plane magnetic field changes the situation. The transformation to the locally 
rotated spin-coordinate frame turns an initially constant external magnetic field, B = B 1, into an inhomogencous 
Zeeman field that accelerates the loss of spin memory of an electron passing through the dot. Below, this effect will 
be described quantitatively. 

For the sake of convenience, we choose the coordinate system (x\, X2) with axes along the crystallographic directions 
ei = [110] and e 2 = [110] such that we can exploit the C 2 „ lattice symmetry of the (001) plane of GaAs. Then, the 
single-particle Hamiltonian in the presence of an in-plane magnetic field B = B 1 with 1m = 1%, 0) can be rewritten 
as 



1 

2m 



Pi 



<J2 

2Ai 



Pi 



2X 2 



(35) 



where Ai,2 = l/(a BS ± a CI ) characterize the length scales associated with the strength of the SO coupling for electrons 
moving along the principal crystallographic directions. Here, &\ 2,3 are Pauli matrices with er 2 = —&2 an d 01 3 = ^iV 
e z = 9^bB\\ is the Zeeman energy associated with the parallel magnetic field. Furthermore, p = P — A is the kinetic 
momentum with P = — iV and the vector potential A = B z [rx e z ]/2 corresponding to an additional perpendicular 
magnetic field. 

In the following, we assume the dot to be sufficiently small to fulfill the conditions £1,2, <C Ai ; 2 and 7, ez 2?Thj 
where £/rh = hD/L 2 is the Thouless energy. 



B. Effect of spin-orbit coupling on weak localization and conductance fluctuations 

For a spin-i particle in a quantum dot connected to adiabatic ballistic contacts, the WL corrections to the conduc- 
tance can be related to the lowest-lying modes of Cooperons in the singlet and triplet channels. The expressions for 
the different Cooperon modes can summarized in the equation 



L ~ 2 



-(tr a L U2G\{e)(J20MGA(e-uj) (36) 



where L, M — 0, . . . , 3. Note that the use of the (perturbative) diagrammatic technique in the description of quantum 
dots is justified if the number of channels in the leads is large, ^> 1. 

In the absence of SO coupling and Zeeman splitting, the Cooperon modes C can be separated into four completely 
independent channels, i.e., one singlet channel Cq = C 00 and three triplet channels Cf 1 — C MM [M — 1,2,3). Or, in 
a matrix representation , 

C = CS, where 5 = 8lm, 
and C obeys the conventional diffusion equation. This leads to the familiar result for the WL corrections j27j, 




However, the SO coupling and Zeeman splitting mix the various components pq ] and, thus, split their spectra. This 
changes the conventional diffusion equation into the matrix equation 

nC(x,x') =<5(x-x')<5, (37) 

where 
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A, 



U = -D[(d Xl +2iA 1 )S-—S 2 ) -D[(d X2 +2iA 2 )5+—Si) +~f5 + ie z r) 



(38) 



Here, S K M = 



-l£ 



KLM 



are spin-1 operators (K,L,M = 1,2,3 ), and s is the antisymmetric tensor. As a 4 x 4 



matrix, S also has zero elements when L = or M = 0. The other relevant matrix i) is defined as 

,L ~~= II <>qm + Im $ql 



V 



(where Ik is the Kth component of ly ) indicating that coherence between electrons with opposite polarization is lost 
on the time scale e^ 1 . 

Eq. ( |3"8| ) is supplemented with a boundary condition at the edge of the dot characterized by the normal direction 
nil = (n 1 ,n 2 ), 



nil ■ (V + 2iA)5 - in x \ 1 S 2 + in 2 X 2 1 S 1 



C = 0. 



(39) 



Because of the boundary condition, the lowest mode of the Cooperon cannot be a mere constant solution, such as 
q = 0. To find the lowest Cooperon mode in this problem, one should make such a gauge transformation (that is, a 
unitary rotation of the Cooperon spin components) 

C = uc'ij- 1 , 



that would transform the original boundary condition into nil • VC = 0. In the new spin coordinate system, one can 
now approximate the lowest Cooperon mode for C by a q = solution, i.e., C = const., and evaluate its eigenvalue 
using the standard methods of the Hamiltonian perturbation theory with respect to the terms generated by such a 
rotation in the initial differential operator H. This program can be realized by applying the transformation 



U = exp 



Si 



A, 



exp[-iip s (x)S 3 ] exp[-iy> (x)], 



where function tp = ip s + ip a transforms from the symmetric gauge to a gauge, where the normal component of the 
vector potential at the boundary of the dot vanishes. This eliminates the lowest orders SO coupling terms from the 
boundary condition, and, thus, (in a small dot, L\ 2 -C Ai^) can be followed by a perturbative analysis of extra terms 
in Eq. (^8|) generated by the rotation U. This results in the 0D matrix equation for the Cooperon, 



7<5 + iezfj 



1 



e™S 3 ) + el(S-S'£) + ef|°(S 2 -^ 



(40) 



which depends on six different energy scales to be discussed in the following. (7 and ez have been introduced earlier.) 
In Eq. (|4C|), the two parameters 



4ttB? 



(\M a(3 \ 2 ) = K E Th (2eB z A) 2 



K E Th (A/X1X2) 



possess the same dependence on the shape of the dot and the disorder in the sample. Here, k is a geometry-dependent 
coefficient. Furthermore, the random quantities M a p are the non-diagonal matrix elements of the magnetic moment 
of the electron in the dot. Note that the difference in the third term in brackets of Eq. (|i(]), containing Tp 1 and e s °, 
reflects the addition or subtraction of the Berry and Aharonov-Bohm phases, as was pointed out in Ref. [ p9|j6C| ]. 
The parameter e\ in Eq. ( ^p| ) is the result of a parallel field induced additional Zeeman splitting, 



el 



2A ^ 



A7 A -i 

i,J=l,2 1 3 



where S, 



LiLj 



a/3 Y P a \ 
j 1 



(XfX 



and are the non-diagonal matrix elements of the dipole moment of the electron in the dot. The quantities Hy 
depend on geometry and the disorder and may be estimated as 5 ~ A/£>rh = g 1 - This yields < £z. A similar 
energy scale has been found in other recent publications J6l[ . 
Finally, the parameter 
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(ii/Ai) 2 + {L2/X2) 2 



«e 



introduces the smallest energy scale through which the SO coupling affects the Cooperon propagator. 

In principle, the form of Eq. ( p!0| ) is applicable beyond the diffusive approximation as it follows purely form symmetry 
considerations. Now, the weak localization corrections can be found from Eq. ( |40| ) as 



g wh oc tr 



C 5-S- 



In the absence of a perpendicular magnetic field (t^ 1 
7, e z , D/\\ 2 < £ Th , this yields 



(41) 



0) and using the assumptions made above, namely 



,WL 



e 2 air 
h 4 



7 + e s l 7 + 4 + 2ef. 



7 + 6^ + 



(42) 



where the fact that e 7 ! <C ez has been used. Furthermore, we introduced the notation 



ai r 



WlNr 

(Ni + N r ] 



The formula in Eq. (Ji^) describes the average tendency of the in-plane magnetoresistance of a dot due to the interplay 
between SO coupling and Zeeman splitting. 

It is interesting to note that the application of a Zeeman field (i.e., an in- plane magnetic field ||]) alone does not 
suppress the weak localization corrections completely as long as ez <C i?Th - whereas in the opposite case, e z > £/rh, 
it does. However, the orbital effect of such a strong in-plane field is already sufficient to suppress weak localization 
as discussed in the previous section. 

In the regime ey°/7, e^/7 — > relevant for the experiments on small dots Q, the form of Eq. ( |42] ) is dominated by 
the first two terms, 



1 



1 



l + e\{B)h 



(43) 



This suggests a possible procedure for measuring the ratio A1/A2. By fitting the experimental magnetoresistance 
data to g WL (-B), one can determine the characteristic in-plane magnetic field B at which the weak localization part 
of the dot conductance gets suppressed by the factor of two. For a dot with a strongly anisotropic shape, this 
parameter would depend on the orientation of the in-plane magnetic field. In particular, B should be measured for 
two orientations of the in-plane field: namely fi[no] for 1 = [110] and i3[Ho] for 1 = [110]. Furthermore, one should 
perform a simultaneous measurement of the two characteristic fields £>| 110 ] and S' 1 j ] in a dot produced on the same 
chip by rotating the same lithographic mask by 90°. The anisotropy of the SO coupling is then obtained directly from 
the ratio 



(%io]#[iio]/ s [iio]#[iio]) - (Ai/A 2 ) 4 , 



that is, independently of details of the sample geometry. 

The results for weak localization part of a two-terminal conductance of a dot and the variance of its universal 
fluctuations are summarized for the limiting regimes in the following equation, 



WL 



2-/3 
2/3S 



air 



(S9 2 ) = 



4/3E 



(44) 



Here, the conventional parameter (3 describes time-reversal symmetry of the orbital motion, s is the Kramers' degen- 
eracy parameter, and £ is an additional parameter characterizing the mixing of states with different spins for strong 
Zeeman splitting. In a small dot where e s °/l ^ as well as £j°/7 — > 0, one obtains the following values for 
the different parameters: (3 = 1 indicates that for weak SO coupling the electron spin splitting cannot violate the 
time-reversal symmetry of the orbital motion; Kramers' degeneracy is preserved (s = 2) for ez < 7, but lifted (s = 1) 
for e z > 7; finally, £ = 1 for ej_ < 7 and Y, — 2 for (e z S>) ej_ > 7 at very strong Zeeman splitting. 
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V. IN-PLANE MAGNETORESISTANCE IN SYSTEMS WITH MAGNETIC IMPURITIES 



In addition to a purely potential disorder, metals and semiconductors may contain magnetic impurities |4C| , |39|| . In 
fact, most of real materials certainly do have them to some extent [ |62"||63fl . In this section, we discuss the influence of 
a dilute magnetic contamination on the quantum transport characteristics of disordered conductors. In particular, we 
describe the suppression of weak localization, and its restoration by an in-plane magnetic field due to a polarization 
of the localized magnetic moments - which slows down the decoherence of conducting electrons and produces an 
observable magnetoresistance, g(B\\) = g c u ss + 5g WL (B\\). 



A. Electron spin relaxation and Korringa time for magnetic impurities 

The coupling between a spin-^ electron and the magnetic impurity spin originates from the exchange interaction 
and can be described by the Hamiltonian 

H s = J S(r)-<7, 

where S is the impurity spin and J is the exchange coupling. Furthermore, a is a vector of Pauli matrices. When 
scattering on a magnetic impurity, the electron can change its spin state. Thus, electrons diffusing in a disordered 
conductor containing magnetic impurities loose their spin memory after the time scale r s determined by the spin-flip 
relaxation rate 

r" 1 ~ v J 2 S{S+ 1) ~ vn s j 2 S{S + 1). 

Here, j the exchange coupling due to a single impurity and n s is the density of magnetic impurities. For dilute 
magnetic impurities, r s ^> r. 

To incorporate spin-flip scattering into the quantitative Cooperon/diffuson analysis of weak localization and UCFs, 
one has to make assumptions concerning correlation properties of this additional source of disorder. Similarly to the 
non-magnetic impurities, we assume the magnetic disorder to be Gaussian (5-correlated with zero mean and variance 

(JS a (r,t)JSp(r',t)) = —^—S af3 6(r-v') 

07W7" S 

at coinciding moments of time. 

There is, however, an important difference between potential disorder, which is static for a fixed configuration of 
impurities in the sample, and magnetic scatterers, which change their spin state after each electron-impurity spin-flip 
event. As a result, magnetic disorder has its own dynamics, and a magnetic scatterer forgets about its initial spin 
state with the so-called Korringa relaxation rate p3] 



-1 „-2 



T 



vn e j 2 — , (45) 



where n Q is the electron density. 

Korringa relaxation, as it is known, takes place due to spin-flip scattering at the impurity of any electron with 
energy |e — ep\ ~ T close to the Fermi level. Therefore, the corresponding relaxation rate is temperature-dependent, 
whereas the single-electron spin relaxation rate l/r s is temperature independent. As a result, for temperatures T > T s 
higher than a certain temperature 

n s 

T s ~ — e F , 

n e 

an impurity changes its spin state faster than any free electron whose individual propagation we may follow, i.e., 
1/tk > 1/t"s- By contrast, for T < T s , the single-electron spin relaxation is faster, l/r s > 1/tk, and the exchange 
field of the impurities can be treated as static when analyzing weak localization effects. 

The relaxation of the impurity spin sub-system should be taken into account in the correlation properties of this 
source of randomness. This can be achieved by saying that, for the same magnetic scatterer, 

(S a (r, t)S p (r, 0)) = 5 aP S(S + 1) e~^^. 
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B. Spin-flip scattering, impurity spin dynamics and polarization of impurities in weak localization 



After an electron flips its spin a few times at encountered magnetic impurities, the coherence between waves 
propagating in clockwise and anti-clockwise directions along the same geometrical path is destroyed. In the regime 
of fast Korringa relaxation, 1/tr > 1/t s , this happens since electron waves traveling towards each other along a loop 
encounter the same magnetic scatterer at different moments of time and, therefore, in uncorrelated initial states. 
In the regime of slow Korringa relaxation, 1/tk < 1/t s , decoherence between these two waves accumulates due to a 
non-commutativity of electron spin operators, whose product would appear due to the coupling to spin-flip transitions 
between the same initial and final states of the same set of magnetic impurities visited in the opposite order. 

As a result, in a magnetically contaminated conductor, the time of flight along trajectories forming the enhancement 
of backscattering is limited by the spin-relaxation time, r s . As in the previous section, a quantitative description of 
weak localization corrections to the conductivity, 



fJ 



WL 



e 

2nh 



dt [C (*;0)-3Ci(t;0)], 



(46) 



requires a diagrammatic evaluation of the singlet and triplet Cooperon correlation functions. In Eq. ([46[), the variable 
of integration t is the difference between moments of time when clockwise and anti-clockwise propagating waves pass 
through the same point in space. As shown in Fig. |l2], when deriving the equation of motion for Cooperons, one should 
include magnetic impurity scattering into both, (a) the self-energy of the impurity-averaged electron Green functions, 
i.e., as a correction to the mean free time, and (b) the two-particle correlation functions as additional impurity lines 
in the ladder diagram. Taking into account the Clebsh-Gordan coefficients that appear upon splitting the Cooperon 
into singlet and triplet channels, the equation for the Cooperon takes the form 



d t - DV 2 



+ y (l + cje-W 7 *) ] Cj(t, f; r, r') = S(t - t')6(r - r'), 



(47) 



where J = (spin-singlet) or 1 (spin-triplet), and Co = 1, c\ = —1/3. The time-dependence accounts for the fact 
that, for clockwise and anti-clockwise trajectories, electron waves test the same scatterer at different moments of time. 
If the relevant time scale determining the value of the integral in Eq. ( fl6| ) is large such that \t\ > tk, the states of 
the same impurity seen by waves propagating in opposite directions are uncorrelated. In this regime, the solution of 
Eq. @ is 



C(q,w) 



1 



Dq 2 — iuj + t s 



(48) 



with the same relaxational pole t s , i.e., the single-particle spin relaxation rate, in all spin channels. This regime is 
realized when tk < t s , that is, at high temperatures T > T s . 

By contrast, for tk > t s , the integral in Eq. ( |4^ ) is determined by the behavior of Cooperons at |£| < tk, where the 
decoherence in the triplet channel takes about thrice longer than in the singlet channel, 



Co(q,w) 



1 



Dq 2 



2t s 



Dq 2 



3 Ts 



The latter feature indicates that a single spin-flip does not completely destroy coherence between clock and anti- 
clockwise paths, given that in both propagation scenarios the magnetic impurity state undergoes exactly the same 
transition. 




FIG. 12. Cooperon for \t\ < tk (both contributions) and |t| > tk (only left diagram). (The thin dashed line represent 
scattering from the non-magnetic impurities while the thick dashed lines stand for (spin-flip) scattering from the magnetic 
impurities.) 
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Eq. ( [47| ) can be solved using the following substitution, 

C;(i 1 /;r,r') = e-^C;(t ) t';r ) r')e^, 

with 

/(*) = ~ (* + c ^ si S ni (! " e~ |t|/TK )) . 

Then the 'modified' Cooperon Cj(t,t';r,r') obeys the conventional diffusion equation (d t — DW 2 )Cj(t,t';r,r') = 
S(t — t')^( r — r ')- After substituting the resulting expression for the Cooperon Cj back into Eq. (f46|), we arrive at 



dO 

Qdj2 



t/t s 



exp 



In the 'low-temperature' limit, where r K x — * 0, 



,WL/ -1 



(0-0 = --(iJr.) 1 



3 exp 



W 3 3/ 2 



v^2 



In- 



_ 2 0+^(l-e-^/ TK ) 
3r s 



<2 = 1, 
d = 2. 



(49) 



(50) 



For comparison, in the limit t k x — > 00, the above equation changes into 



WL 



(r K 1 ^oo) = -^-( J Dr s ) 1 - 



2 

ln^ 



d = l, 
d=2. 



The difference between the dynamic and static results in Eqs. (^l|) and (^) manifests itself in a weak temperature 
dependence of the conductivity, 



a(T > T s ) - a(T « T s ) a — (Dr.) 1 "* 



1.2 
1.9 



rf= 1, 
d = 2. 



(51) 



So far in this section, we addressed only the suppression of weak localization by spin-flip scattering, whereas our main 
goal was to identify possible effects of an external magnetic field leading to an observable in-plane magnetoresistance. 
An in-plane magnetic fields tends to polarize the paramagnetic impurities, thus, opening a gap = g( imp ^ hbB\\ 

for the electron spin-flip relaxation process. At high magnetic fields, e% mp ^ > T, the allowed electron energies are 
insufficient to flip the impurity spin, i.e., spin-flip scattering is suppressed. Thus, at B« > B s , where 



B S (T) = T/(^ imp ) 



(52) 



magnetic impurities act as non-magnetic 'potential' scatterers, with the only difference that their potential profiles 
are different for spm-f and spin-| electrons. Treating the spin-| and spin-| electron subsystems separately, one may 
conclude that, for T <C e^ mp \ each of them contributes a 'restored' value towards the weak localization correction to 
the conductivity (conductance), which results in its overall drop across a broad range of in-plane magnetic fields, 

I ^ 

-,2 



ln^, 



2D: 



Q1D wire; 
OD dot. 



(53) 



In Eq. 



r w > r s is a 'true' inelastic decoherence rate due to the electron-electron interaction or external electro- 



magnetic noise. 

The quantitative analysis of a system with fully polarized static magnetic impurities implies the evaluation of all 
Cooperon channels, Co, C\q, and Ci±i- The first two happen to have a gap in the spectrum, both due to i) the splitting 
between up/down-spin bands caused by the combination of B» and the mean exchange splitting due to impurities, 
and ii) the difference between the scattering conditions for up/down-spin electron. By contrast, C\\ and Ci_i remain 
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massless. When t v 3> t s , the crossover between the two limits in Eq. (|5^) can be roughly described by substituting 
the electron spin-flip rate at intermediate fields, 

= -^—^ x where u = 5||.9 (imp V B /T = B\\/B B , (54) 

into Eq. (|50|). The field-dependent electron spin- flip rate in the latter expression manifests the activational character 
of the impurity spin flip process and was calculated after taking into account thermal occupancy of spin-split impurity 
states. For a2D electron gas or a thin metallic film, this yields a magnetoconductivity with the characteristic form, 

e 2 fifV) 2 B n <B, 



ff(S||)-ff(0) = T ln-— «- T x< "V*; " (55) 



^ 1 §7 B ll>^ 



C. Effect of impurity spin dynamics on mesoscopic conductance fluctuations 

At B\\ = 0, contamination by magnetic impurities also suppresses mesoscopic conductance fluctuations - as a result 
of the combination of electron spin-flip scattering with the impurity spin dynamics. In a small sample, the mesoscopic 
part of the interference correction to the conductance, Sg(t), is finite for any instant configuration of magnetic 
scatterers, but it is specific for each particular snapshot of the magnetic subsystem. As the localized magnetic 
moments undergo a random temporal evolution, this interference correction to the conductance also fluctuates in 
time, thus, leading to a self-averaging of the UCFs to much smaller values. By contrast, in the regime of strong 
magnetic fields, Bu 3> B s , the impurity spins stay polarized, such that the random potentials for spin-f and spin-| 
electrons are static (though different), and UCFs are fully restored. 

To describe the reproducible part of conductance fluctuations in a DC-current measurement over the entire range 
of fields, one has to analyze the time-averaged conductance, 



and its variance, 



I f T 

9 = lim ~ / dt g(t) = (g) + Sg, 



(Sg 2 ) = T Bmi/ dt (Sg(t)Sg(O)) = lim (6 g(t) 5 g(0)). 



q L f '^<J 

The correlation function (Sg(t)Sg(0)) determining (Sg 2 ) can be expressed in terms of diffusons T>j(i~i,r]';t;r,r') with 
an additional time variable t and Cooperons as in Eq. ( ]47| ) in the usual way - via the standard set of perturbation 
theory diagrams. 

The diffusons T>j(ri,ri';t;r,r') describe the quantum diffusion of electrons across two distinctive time intervals 
separated by t. They obey the following equation, 



I (l - cje-KI/^) | Vj{ V , 77'; t; r, r') = S( V - V ')S(r - r'), 



fiL - DV 2 + - 1 - c,e 

T s \ 

where the loss of spin memory of the magnetic impurities due to the Korringa relaxation is taken into account. In 
the relevant limit of |£|/tk — > oo, where e 1*1/ — + 0, the same simplification occur as the ones leading to Eq. ( |4g| ) 
for the Cooperon: the poles of all (singlet and triplet) diffuson components acquire the same relaxational gap equal 
to the single-electron spin- relaxation rate, i.e., 

X>(q,w)= ' 



Dq 2 — iuj + t s 1 



As a result, one arrives at strongly suppressed fluctuations. For a wire or dot geometry, their variance for a given 
value of the magnetic field follows the field-dependence of the spin-flip scattering rate for an electron in Eq. (M) , as 



(Sg 2 



2 



X 



(f s /7-fi) 3/2 wire geometry, 
(t s /t csc ) 2 dot geometry, 
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where t s = t s (Bh), t osc is the electron escape rate from a dot defined in the previous chapter, and tr = L 2 /D. 

The suppression of mesoscopic fluctuation in magnetically contaminated conductors as well as their restoration 
by a magnetic field have been observed, both in metallic ^| and semiconductor |40| microstructures. It is worth 
mentioning that this characteristic behavior of weak localization and mesoscopic fluctuations may be used for testing 
the nature of decoherence |35|,|36| even in purified materials - in particular, when its suspected origin is spin-flip 
scattering at residual magnetic impurities p3,K3[. 



VI. 2D LAYER SYMMETRY, IN-PLANE MAGNETORESISTANCE AND THE NON-LINEAR 

SIGMA-MODEL 



In the previous chapters, we have shown how to analyze "simple" problems using the diagrammatic perturbation 
theory technique. In this section, we give an example of a problem that is more convenient to tackle with the 
supersymmetric (SUSY) sigma-model method: the analysis of the interplay between spatial symmetry of a quasi-2D 
electron gas with few occupied subbands and an in-plane magnetic field.^ The main results have been discussed in 
chapter [II. Here, we first give a brief overview of the SUSY technique and, then, present the derivation of the SUSY 
sigma-model action for a multi-subband 2D system subject to an in-plane magnetic field. Finally, we discuss the 
properties of the Cooperon matrix and the consequent implications for the magnetoresistance. 



A. The field theoretic technique 

In Sec. |n] the diagrammatic approach to calculating correlation functions in the perturbative regime has been 
introduced. An alternative approach is provided by the coherent state path integral. The (retarded) Green function 
can be represented as a field integral, 

G+(r,r';e) = (r|(e + ^0 - #) - V) = ^ J DsDs* s*(r)s(r') e l I drs *<- £+ 

where Z — J DsDs* exp[i J dr s* (e + —7i)s], and H = p 2 /(2m) + U(r). Here V represents the impurity potential 
which is assumed to be drawn from a Gaussian white noise distribution, cf. ([}]). 

Unfortunately in this form it is not possible to carry out the disorder averaging: Due to the presence of the 
partition function as a normalization factor, Z _1 [V], the random potential appears in the numerator as well as in the 
denominator. One method to circumvent this problem is supersymmetry p7{ . When one is considering only single- 
particle properties of a system, there are two equivalent formulations of the path integral, namely by using bosonic 
or fermionic fields. Supersymmetry now exploits the following property of commuting (s) versus anti-commuting or 
Grassmann (x) variables: 

J ds* ds e~ a ' Ma = det -1 M, J d X d X e~ mx = det M. 

Thus, combining both variables into a 'supervector', ip T = (s, x), yields the result J dip' dip exp[— ip' M (g> l bf ip] = 1, 
where the superscript 'bf ' stands for 'boson-fermion'. Applying this to the partition sum, it is automatically normalized 
to unity, 2 = 1, and the impurity averaging is straightforward. The evaluation of a two-particle correlation function 
requires the introduction of two sets of fields, covering the advanced and retarded sector. With ip T = (si, xi, s%, X2) 
the correlator can be written as 

(G + (6+|)G-( e -|)) =-Q D^s^s^e-^^-^-^y 

where lu + = u + iO and a^ R is a Pauli matrix in advanced/retarded space. Furthermore, ip = ip^ L with L = 
<73 R (g>-Ebb + l aT ®Eff , where -Ebb and Eg are projectors onto the boson-boson and fcrmion-fermion block, respectively. 



4 In fact, the present, perturbative problem is accessible by diagrammatic methods as well. However, to get the interplay be- 
tween inter-band correlations and disorder scattering reliably under control, the formalism of field integration has the advantage 
that the fully microscopic aspects of the problem are processed in the early stages of the derivation |67[. 
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By introducing a source term, Sj = — J dr (J^ip + V>J), different correlators of Green functions can be obtained from 
the generating functional 2[J] by taking derivatives with respect to the source field J. In the following, we will 
suppress the sources and consider only 2[0]. 

Now the impurity averaging of the partition function leads to a quartic term in the fields ip, 



/ e i/dr^\ =exp r__L_ f drMf 



By Fourier transformation to momentum representation, one can identify the slow modes, 

dr{^) 2 = (V'piV'pJfeaV'pJ 

(W'pV'-p+qX^-p'V'p'-q) + C^pV'-p'XVV-qV'-p+q) + C^p'-qX^-p+qV'-p')) > 



p>p';q 
where |q| <l _1 . 

The first term corresponds to slow fluctuations of the energy which can be absorbed by a local redefinition of the 
chemical potential. Thus, we concentrate on the remaining terms: The second term generates the diffuson contribution 
while the third term yields the Cooperon contribution. Enlarging the field space^Jby defining ^ T = (ip T , i/^)/\/2, the 
last two terms can be rewritten into a single contribution, Str C( c l)C( — where £(q) = J2 P ^(P - q)^ 1 ! - P)- The 
components of the newly defined vector fulfill the symmetry relation ^ = (C^) 7 , where C — a^^E^ + ia™ ®Efj. 
This symmetry corresponds to time-reversal (ip — > ip* , Ti. — ► 7i T ). In the absence of the symmetry-breaking energy 
difference, u> = 0, the action is invariant under rotations — > U^> , where ULU^ = L and U T — CWC T . Thus, 
U € Osp(4|4). 

As a next step the quartic interaction is decoupled by a Hubbard-Stratonovich transformation, introducing the new 
(supermatrix-) fields Q: 



exp 



1 



AlTVT 



dr 



DQ exp 



^ jdr StrQ 2 -^ l<h-mW 



where Str M = tr A/bb — trMgf. The symmetries of Q reflect the symmetries of the dyadic product \& <g) namely 
Q = CLQ T (CL) T . Now the resulting exponent is only quadratic in the original 'J'-fields. Therefore, the Gaussian 
integral can be readily evaluated, yielding the action 

S[Q] = J dr Stl 'Q 2 + 2 J dr Str ln ^ _1 ' (56) 

1 A2 , i 



where (Z) = fDQ exp(-5[Q]) and Q~ x = ^p 2 - e F 



- c F -r — <J 3 t 27 



To extract an effective low-energy, long-wavelength field theory from this action, a saddle point analysis has to be 
performed. Variation of (|5^ ) with respect to Q yields Q sp (r) = iG(r, r)/(-nv). Neglecting the small energy w, the 
Ansatz Q sp constant and diagonal leads to 



Qsp = -- I dti - r^- = sgn(Q sp ). (57) 



Thus, the saddle point Q sp has the meaning of a self-energy. Analytic properties of the Green function single out the 
solution <5s P = C3 11 . 

In fact, the action is invariant under transformations Q — ► TQT^ 1 , where T constant: instead of one saddle point 
one obtains - at u> = and in the absence of symmetry breaking sources - a degenerate saddle point manifold Q 2 = 1. 
Fluctuations around the saddle point can be subdivided into longitudinal modes, [SQi, Q sp ] — 0, and transverse modes, 
{SQ t ,Q sp } = 0. The longitudinal modes 5Qi leave the saddle point manifold Q 2 = 1. Therefore, they are massive 
and do not contribute to the low-energy physics of the system. In the following, we concentrate on the transverse 



5 Generally, each discrete symmetry leads to a doubling of the low-lying modes and, thus, should be incorporated by doubling 
the field space p8|. 
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modes 5Qt- The parameter which stabilizes this distinction is k^t ~ epr, i.e., the following considerations are valid 
in the quasi-classical limit. We proceed by expanding the action around the saddle point in the slowly varying fields 
Q(x) = T(x)Q sp T~ 1 (x). Separating the fast and slow degrees of freedom with p — > p + q, this expansion yields 



S ~ i / dr I dp Si ! 



1 



(SoT-ipqT) 



2 " J " 2m 2 

The integral over fast momenta, p, can be performed using the following representation for the Green function, 



2 



i^(l + S( T 3 AR )Gg( P ) 



Then, J dp £?o(p) = — i7w<73 R , and 

/" dp Str [0o(p)T- 1 p4T0o(p)T- 1 p4T] = if / ^ G o ^ Str [(1+saf jT" 1 qT(l-aof qT] . 

V v ' 

= m 2 irvD/2 

Finally, summing over s and using the cyclic invariance under the trace, the effective action takes the form of a 
non-linear a model, 



S[Q]=~y J dr Str [D{dQf + 2iuj + cr : 



(58) 



Note that the effect of a weak magnetic field is to generalize the derivatives to d — d — iA[a^ R , .]. For the quasi- 
twodimensional system in an in-plane magnetic field this will be discussed in more detail in the following. 

In the perturbative regime, an expansion of the effective action around the saddle point in the generators W, where 
Q = e~ w ^ 2 a^ R e w ^ 2 and {W, cr^ R } = 0, reproduces the diagrammatic results. By contrast, in the non-perturbative 
regime, the action is dominated by zero-modes which require an integration over the whole saddle point manifold 
Osp(4|4)/(Osp(2|2)®Osp(2|2)). For our purposes, a perturbative expansion will be sufficient. 



B. Magnetoresistance in a multi-subband electron layer with a possible Berry-Robnik symmetry 

As pointed out earlier, the magnetoresistance in quasi-twodimensional electron systems subject to parallel fields is 
very sensitive to the presence or absence of T^-symmetry. While in chapter ^ we restricted ourselves to a qualitative 
discussion of the effect, here the quantitative results are presented. To explore such type of phenomena one needs 
to construct an approach which on the one hand is sensitive to microscopic details in z-direction while on the other 
hand should be capable of efficiently describing large scale in-plane properties. This task can efficiently be addressed 
within a field integral formalism. The starting point of the derivation is a supersymmetric field integral with action 

Sty] - i J d\^(e F - w+af + ±-{dl + {d v -iBzf + d 2 z ) - W(z) - V(x,yj)i/>, (59) 

where W is the confining potential of the 2DEG, and V a disorder potential. To simplify the analysis, it is assumed that 
the disorder potential does not depend on the z -coordinate. Given the typical architecture of 2DEGs, this is certainly 
a justified zeroth order assumption. In, e.g., high mobility 2DEGs (or 2DHGs) in GaAs/AlGaAs heterostructures, 
the mobility is limited by a long-range random potential, created by charged impurities located far from the plane. 
Later on, this condition will be relaxed by generalizing V according to V(x,y) — » V(x,y) + U(x,y,z), where U is 
weak and can be treated perturbatively in a sense to be specified below. 

The confinement in z-direction is responsible for the size quantization which entails a subband structure of the 
system. To make progress with the action, Eq. (|9|), an orthonormalized set of wavefunctions {4>k\, diagonalizing the 
z-dependent part of the problem, is introduced: 

(-^-9 2 + W(z)- e k ) fe = 0. (60) 

Expanding the original fields ip in the complete set of eigenfunctions 4>k, that is ip(x,y,z) = ^2j,ipk{x,y)(j)k{z) 1 the 
action takes the form 
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Sty] = i Jd 2 r 4> k ( [e F - w+of - e k + || - y (as, y)] fc fc , + ^ ([d y - iAf) ^ ) fa , 
where the integration extends over the x-y-plane, summation over k, k' is implied, and 

A k k'=B dzcf) k (z)zcf)k'(z) (61) 



is the vector potential. I.e. the magnetic field couples to the dipole matrix elements d k k' = Jdz (f>k(z)zcf>k l (z) that 
contain detailed information about the microscopic symmetry properties of the system. 

Using the same steps as explained in the previous section , the action in terms of the slow supermatrix fields Q kk > 
reads 

S[Q] =~J d 2 r StrQ 2 + ± J d 2 r Str In L F -w + af-e+±{<% + {d v -iAal*Y)-^C^ , (62) 

where a compact fc-index free notation has been introduced. Here the energy matrix e = diag(eo, ei, . . .) contains the 
subband energies. The next step in the construction of the effective theory is the saddle point analysis. Functional 
differentiation of the action with respect to Q obtains the equation 

A k = — { d 2 p ' 



m>J iSar-eF-ek + ^ + ^Ak 



for the diagonal elements Ak of the saddle point matrix Qkk' ■ At this stage one has to specify the relative position 
of the Fermi energy E-p and the subband energies ek- Below, we will explore the case where M bands with energy 
e k < E-p (k = 0, . . . , M — 1) exist. This leads to A k — a^ R for k < M, and A k = otherwise, where it has been 
assumed that the highest occupied subband cm-i lies well below (farther than r _1 ) the Fermi level. Based on this 
solution, the low-lying fields of the theory can be represented as Q — TAT -1 , where A = {A k 8 kk i}, and the final 
expression for the general slow action reads ]15|,[16[ 



S[Q] = ~Y j dS ^ Str { 4tLja 3 R Qk + D k (d k Q k ) 2 ) + ^-JdS Str (Xkk-vTQ k a7Qk), (63) 

h ' k,k' 

where dk = d — ie y Akk [^3 R , ■ ] ■ The sum ^2 k involves only the occupied subbands k = 0, . . . , M — 1. Furthermore, 

X kk > = \{D k + D k ')-, \^—A kk ,A k , k (l - 5 kk >). (64) 

2 (e kk 'Ty + 1 

Here D k is the diffusion constant of subband k, and e kk > = e k — e k > . The first line of ( p3| ) gives the conventional result 
of a 2d system while the second line describes the coupling of the subbands induced by the magnetic field. 

To prepare the one Cooperon approximation to the conductivity, one has to expand the fields Q to lowest non-trivial 
order in some generators. It is convenient to decompose the generators into 'diffuson' (d) and 'Cooperon' (c) blocks. 
It is clear from the structure of the second order action, that it does not couple between the l d' and the 'c' sector, i.e., 
S = S d + S c . The diffuson action S d is generated by those fields that commute with ct™ and, thus, do not couple to 
the magnetic field, while the Cooperon action S° is generated by those fields that anticommute with tr™ and, thus, 
are field-sensitive. 

Here we are only interested in the Cooperon action being responsible for WL corrections. The kernel appearing 
in S c is the 'inverse of the Cooperon'. More explicitly, the Cooperon C, which in our formulation is a matrix in the 
discrete space of fc-indices and diagonal in q-space, is obtained by inverting the matrix 

(2zlu 2 i 2 

--p— + (q-2A fcfc ) 2 + — / Xkk" Skk' H — 7= = 
E>k k / \JDkUk' 

The magnetoconductance is determined by the specific form of this matrix. The conductivity is given as a = a + Act 
with ct° = X^fcLo 1 a ki wriere o'fc = vDk is the Drude conductivity of subband fc, and 

2 M_1 

A ( r = -- {Cw= )kk- (65) 

k=0 q 
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In general, the field-dependent terms will render C massive, i.e., the weak localization corrections will suffer from a 
field induced suppression. However, there is the situation mentioned above, where z-inversion, V z : z i— > — z, is an 
(approximate) symmetry of the Hamiltonian, \H, V Z \ ~ 0. 

For systems with an exact 'P 2 -symmetry, the eigenfunctions <t> k obey V z (j>k = {~) k <t>k- The definition of the vector 
potential matrix A then implies 



A 



kk 1 



Akk' k + k' odd, 
k + k' even, 



and, thus, the same holds true for X kk i ■ 

This structure bears consequences on the Cooperon mass. To analyze this point, consider the spatial Cooperon 
zero-mode, 



2 (^JJ~ ^2 Xkk"Skk' 



it o 1-0 )/.•/••< = 2 ! TT zL •V/.A-W- + ~^=^=\ ■ (66) 



This matrix has determinant zero implying that there is a Cooperon mode which is not affected by the field. In fact, 
it is straightforward to verify that the M-dimensional vector 

X = M- 1 / 2 J2(-) k VDk~e k , (67) 

k 

where J\f = ^2 k Dk, is annihilated by Cq As the zero-mode matrix is symmetric, one can, in principle, construct 
a complete set of orthonormal eigenvectors, {Xo = X, Xi, . . . Xm-i}, with eigenvalues {0, Ai, . . . , Aa/_i}. 

Furthermore, due to Akk — 0, the full Cooperon kernel is separable (i.e., it is the sum of a spatial and an 'internal' 
operator). Inserting the result into Eq. ( |65| ) yields 

Aa(B) = -~J2^—. (68) 

A:;q 

This is our final result for the conductivity. Notice that, due to Ao = 0, the weak localization corrections do survive 
the magnetic field; carrying out the q-summation leads to the usual logarithmic correction to the Drude conductance. 



Thus, even at high magnetic fields, a logarithmic temperature dependence - see Chap. Ill - of the conductance should 
be observable. All other eigenvalues are proportional B 2 and, thus, display the usual field dependence. 

In the following, let us concentrate on the behavior of the lowest eigenvalue if the system is not exactly inversion 
symmetric. An asymmetry can be caused cither by the confining potential or by a z-dependence of the random 
impurity potential. 

In the case of an asymmetric confining potential, Eq. ( |66| ) generalizes to Akk 1 — * Akk' + 5 Akk' , where 8 Akk' is 
assumed to be much weaker than the symmetry allowed elements Akk' , k + k' odd. Similarly, there are non- vanishing 
but small matrix elements 8Xkk' for both k + k' even and odd. To lowest order in perturbation theory, the presence of 
these matrix elements shifts the zero-mode eigenvalue Ao (q) of the unperturbed Cooperon mode at momentum q by 
the amount <5Aq (q) = X T 5Cq : X, where 1 is the perturbation contribution to the Cooperon operator. Explicitly, 

{5C^) kk , = [-Aq y SA kk +ASA 2 kk +2j2^kk"]Skk'+26X kk ,. (69) 

k" 

Combining these equations and making use of the definition of the zero-mode eigenvectors d67j ) yields 

SX ( o S) (<7mi„) = DkDk ' ( SAkk - 5Ak ' k 'f + Jf SXkk ' ■ 

k,k' k+k' even 

Thus, the Cooperon acquires a mass term ~ B 2 . 

The influence of z-dependent impurities has a similar effect. A potential with a generic z-dependence will not be 
inversion symmetric, implying that, somehow, the Cooperon must pick up a mass. Assuming a Gaussian distributed 
potential U(r, z), 

(C/(r, z)) - 0, (U(v, z)U(r', z')) = 7 2 <5(r - r')5(z - z'), 
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for small 7, it is sufficient to consider the lowest order non-vanishing contribution in U to the action. We arrive at 



S[Q] = S [Q] + (^ff f dS^Tkk, Str (Q fe (r)Q fc ,(r)), (70) 

^ kk' 

where the coefficient Tkk'k"k>" = j dz 4>k(z)<fik' i z )^k" (z)<f>k"' ( z ) > an d ^kk' = ^kkk'k' is positive. This expression tells 
that the z-dependent scattering tends to lock the fields Qk- For 7 large, only field configurations {Qk = Q} with 
no fc-dependence survive. The physical mechanism is the following: Scattering in z-direction leads to a coupling 
between the different k -bands. Thus, the formerly independent diffusons and Cooperons are coupled, too. The 
formerly massless Cooperon channel does not survive this coupling as the A;-space eigenvector X associated with 
the eigenvalue A = is staggered in k, cf. Eq. (|67|) , i.e., it stands orthogonal on the field configurations that are 
compatible with the locking. 

If the coupling due to the impurity scattering is smaller than the field induced subband coupling, the shift of the 
lowest eigenvalue is again obtained by first order perturbation theory. Then, 

st p) = ^ E ( 71 ) 

fc+fc' odd 

Or, (5A imp) - 1 /(N't'), where r' has the meaning of a scattering time perpendicular to the plane, i.e., between the 
subbands. 

This result which does not depend on the magnetic field holds true only for sufficiently large fields. For smaller 
fields, the disorder induced mass term fixes the preferred eigenvector. To compute the mass of the completely locked 
Cooperon, consider A; = X^C^_ X/, where the Cooperon operator is given by Eq. (|66|), and the 'locked' vector X; 
reads X/ = A/" -1 / 2 Y^k y/Dk &k- Explicitly computing the matrix element leads to 

A; = 772 ZJ D * D * ( Akk - Ak>k ' ) 2 + jf XJ Xkk> ■ ( 72 ) 

At low fields, the mass of the Cooperon increases quadratically with B according to Eq. J7^), but then, due to Eq. j7l|), 
it levels off at large fields. The characteristic field B c can be estimated by comparing Eqs. ( |7l| ) and ( |72| ) which yields 
B c ~ E/vy^Jt/t' /d, where E stands for the typical energy separation between subbands and d sets the scale for the 
width of the quantum well. 

As the expressions obtained above are rather lengthy, it is helpful to consider some specific examples. We concentrate 
on the experimentally most relevant case M = 2 and, for simplicity, choose^ Dq = D\ = D. Diagonalization of the 
2x2 Cooperon matrix, 



n-X _ ( (q — A) 2 + JjXq\ JjX l 

~\ %X m (q + A) 2 + 



yields 



X = q 2 + A 2 + ^X 01 ±2^(Ac i r + ^Xl, (73) 



where A = Aoo — An, and Xqi — DA^ 1 /(1 + (eior) 2 ) obtains from (|64|). The corresponding magnetic decoherence 
times read 1/rg = -DA. Note that at small magnetic fields, 2Aoi <C I/t - ^, the symmetry mechanism is ineffective. 
Irrespective of A, the magnetoconductance yields 

e 2 

a(B) - (j(0) ~ 2— XojTj, 

which shows the usual low-field quadratic dependence on B. However, the coefficient is diminished by the factor 
f/(f + ( £l0 r) 2 ). 



Admitting for different diffusion constants Do 7^ Di does not change the results qualitatively. 
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At large magnetic fields, 2Ani 3> if the confining potential is fully symmetric (A = 0), the result reduces 

to 1/rg = 2Aoi(l ± f). While 1/tb = 4Ani leads to a logarithmic field dependence (see Eq. (|3C|)), due to the field- 
insensitive 1/tb = 0, the conductance maintains its temperature dependence through even at large fields. A slight 
asymmetry of the confining potential entails a finite A, which leads to l/r^iB) ~ DA 2 + 2X 01 (1 ± f). Thus, the 
temperature dependence remains as long as DA 2 < 



O Aa(T) saturates 

• Aa(T)~lnT 

□ Aa(B)~B 2 

M Aa(B)~ In B 

asymmetry 

FIG. 13. Different regimes of B- and T-dependence of the weak localization corrections. 
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The regimes with different field and temperature dependences of the conductance are shown schematically in 
Fig. [l| p^Jl^| . The fields B^ and B^ are defined through 1X m T$ = 1 and DA 2 t$ = 1, respectively. I.e. 

B = 1 B * = 1 

where dkk' are the dipole matrix elements defined above. 



VII. SUMMARY 



In summary, we presented three mechanisms of quantum in-plane magnetoresistance in two-dimensional electron 
systems, such as quantum wells, heterostructures, and inversion layers in field transistors, in lateral quantum dots 
prepared of these materials, and in thin metallic films. These mechanisms include: the purely orbital effects possible 
due to subband mixing in a 2D structure by the magnetic field, the effect of an interplay between the spin-orbit 
coupling in the 2D electron dispersion and Zeeman splitting by a magnetic field, and the field effect on the efficiency 
of spin-flip processes in materials contaminated by paramagnetic impurities. 
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